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PROCESS, DEVICE AND SOFTWARE PACKAGE FOR SMOQTHJN^ A SUBSURFACE 



The invention relates to a process for smoothing a 
subsurface property in a geological structure represented by 
seismic measurements. 

The invention also relates to a device for smoothing a 
subsurface property in a geological structure represented by 
seismic measurements. 

Finally, the invention relates to a computer software 
package permitting the operation of a programmable device for 
smoothing a subsurface property in a geological structure 
represented by seismic measurements. 

The document WO 01/63323 Al discloses a process for 
processing seismic data comprising the steps consisting in: 
obtaining a volume of seismic data covering a predetermined 
volume of ground; determining for each voxel of the volume of 
seismic data the local orientation of the seismic data; 
determining for each voxel if there exists an edge in the 
neighborhood, and executing a smoothing operation on each voxel 
in the volume of seismic data, in which the direction of the 
smoothing operation is the local orientation of the data, and in 
which the operation of smoothing does not exceed the edge, so as 
to obtain a volume of processed seismic data, in which the 
quantities associated with each voxel in the volume of processed 
data is the result obtained by the execution of the operation of 
smoothing in the voxel in the volume of seismic data. 

The document WO 02/13139 Al discloses a process for 
processing seismic images comprising the steps: of obtaining a 
set of data of an initial bi-dimensional or tr i-dimensional 
image, in which element of the set of data is the initial image 
intensity of the point of the image; computing for each point 



PROPERTY 



the partial derivatives of the element in the n direction to 
obtain a set of values derived from the partial derivatives; 
computing for each point a symmetrical structural square matrix 
from the values of the partial derivatives; carrying out a point 
5 by point iteration weighted by a variable near zero when the 
point is near an edge and near 1 when it is very far from an 
edge; and repeating these steps a number of times to obtain the 
processed image. 

A first object of the invention is to improve the smoothing 

10 of the subsurface properties in a geological structure 
represented by seismic measurements. 

A second object of the invention is to permit a simple and 
rapid smoothing of the subsurface properties, whilst not 
smoothing the discontinuities. 

15 The invention has for its object a process for smoothing a 

subsurface property in a geological structure represented by 
seismic measurements, in which there is constructed a continuous 
function S ijfk (t) by interpolation or approximation of the 
discrete seismic traces of a multi-dimensional seismic matrix, 

20 said function being designated by "continuous local seismic 
trace", comprising the following steps: 

a) . using as optimum offset hij, pq , k of two adjacent 
continuous local seismic traces S ij/k (t) and S pq , k (t), the offset 
value rendering maximum their correlation function; 

25 b) taking as the conditional neighborhood of a central 

continuous local seismic trace S ij/k (t) the sub-neighborhood 
consisting in adjacent traces S pq , k (t) corresponding to optimum 
offsets hij /Pq , k associated with correlations R±j, pq , k (h ij/Pq , k ) 
greater than a predetermined threshold comprised between 0 and 

30 1; 
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c) . choosing a subsurface property to be smoothed in the 
conditional neighborhood of a point (i, j, k) of a reference 
"central" continuous local seismic trace; 

d) . offsetting the subsurface properties of the conditional 
neighborhood by translating a current variable of the value of 
optimal offset hij, pq , k ; 

e) . taking as the smooth value at the point (i, j, k) an 
average of the offset subsurface properties from step d) . 

According to other characteristics of the invention: 

- the average of the surface properties of step e) is an 
average which can be weighted, for example by the maximum 
correlation value corresponding to the optimum offset, 

- the average of the subsurface properties of step e) is 
selected from the following set: arithmetic average, geometric 
average, harmonic average, 

- there is chosen as the subsurface property to be smoothed 
the amplitude reflected and detected by geophones, 

- if necessary, the set of smoothing processes is repeated 
a certain number of times to include the smoothing, 

the multi-dimensional matrix cross-sections of smooth 
properties are visualized on a visualization screen. 

The invention also has for its object a device for 
practicing the process according to the invention, comprising 
means to use as optimum offset of two adjacent continuous local 
seismic traces the offset value rendering maximal their 
correlation function, means to offset the subsurface properties 
of the conditional neighborhood by translating the current 
variable of the value of the optimum offset hij, pq , k , means to 
select a subsurface property to be smoothed at the conditional 
neighborhood of a point (i, j, k) of a reference "central" 
continuous local seismic trace, means to offset the subsurface 
properties of the conditional neighborhood by translating the 



current variable of the value of the optimum offset h ij#pq , k/ and 
means to take as the smoothed value at the point (i, j, k) an 
average of the subsurface properties offset in step d) . 

According to other characteristics of the invention: 
5 - the device comprises memorization means and means for 

visualizing seismic parameters determined with the help of the 
process according to the invention. 

Finally, the invention has for its object a computer 
software package, comprising program code elements to execute 
10 the steps of a process according to the invention, when said 
program is executed by a computer. 

The invention will be better understood from the 
description which follows, given by way of non-limiting example, 
with reference to the accompanying drawings, in which: 
15 - Figure 1 represents schematically a three-dimensional 

seismic matrix and a continuous local seismic trace. 

- Figure 2 represents schematically an example of local 
neighborhood Nij /k consisting of a set of continuous local seismic 
traces S pq , k (t) adjacent the central reference trace Sj_j /k (t) 

20 itself located on the vertical of the geophone Gij and 
interpolating the seismic amplitudes about t=k. 

- Figure 3 represents schematically a vertical section of 
conditional local neighborhood N ijfk (r) of a reference continuous 
local seismic trace Sij, k (t). 

25 - Figure 4 represents schematically a vertical cross- 

section of a three-dimensional seismic matrix with an optimum 
offset hij, pq , k and a conditional neighborhood Nij #k (r). 

Figure 5 represents schematically a vertical cross- 
section of a three-dimensional seismic matrix analogous to that 

30 of Figure 4, with the presence of a fault or discontinuity. 

- Figure 6 represents schematically a functional organigram 
of a process according to the invention. 



With reference to Figures 1 to 5, identical or functionally 
equivalent elements are designated by the same references. 

In Figure 1, a three-dimensional seismic matrix is obtained 
by taking measurements recorded by geophones Gij disposed on a 
5 network x,y at coordinate points i, j. The taking of these 
sampled measurements according to time is represented on a 
descending axis t representative of the depth or of the vertical 
descent from the surface of the ground or the sea. The 
measurements are characterized by their amplitude, for example 
10 an amplitude taken by the geophone Gij at the time or at the 
depth of sampling t k . The discrete measurement carried out by 
the geophone Gij in time or at depth t k is called the seismic 
amplitude Sij, k . 

The set of seismic amplitudes corresponding to a geophone 
15 Gij of coordinates i, j is a uni-dimensional matrix (Siji, 
Sij 2 ,...., Sijic, - - . . , S ijN ) called a discrete seismic trace, 
because this uni-directional matrix corresponds to the trace 
according to the point of horizontal coordinates i, j of the 
three-dimensional seismic matrix obtained by seismic 
20 measurements. 

The vertical axis t oriented along a descending vertical 
line usually designates time, but can also be considered to 
represent a depth from the surface. 

The invention also relates to the application of a third 
25 coordinate t representative of time, as well as a third 
coordinate t representative of depth. 

From the discrete seismic trace located on the vertical of 
a geophone Gij, there is defined by interpolation or 
approximation discrete values about t=t k =k, a continuous function 
30 Sij, k (t) which is designated as "continuous local seismic trace". 
The methods of approximation or interpolation of discrete values 
to give rise to a continuous function are numerous, and comprise 
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particularly polynomial interpolations or approximations, as 
well as polynomial trigonometric interpolations or 
approximations . 

Any other variant of interpolation or approximation within 
5 a continuous function can also be used in the present invention 
to provide a "continuous local seismic trace". 

In Figure 2, several continuous local seismic traces define 
an example of neighborhood of a reference "central" continuous 
local seismic trace Sij, k (t). The neighborhood of a continuous 
10 local seismic trace Sij /k (t) is defined as the set of continuous 
local seismic traces whose horizontal spatial indices p, q are 
adjacent the horizontal spatial indices i, j of the reference 
continuous local seismic trace. 

By way of example, the horizontal spatial coordinates p,q 
15 corresponding to the geophone G pq are adjacent the horizontal 
spatial coordinates i, j corresponding to the geophone G±j if the 
absolute values of the differences i-p and j-q are less than 
given whole numbers, for example 2. 

In this case, as shown in Figure 2, the continuous local 
20 seismic trace Sij, k (t) is associated with eight adjacent 
continuous local seismic traces surrounding the "central" 
continuous local seismic trace Sij, k (t). 

In the case of continuous local seismic traces produced by 
seismic measurements, the profile of the geological horizons 
25 introduces vertical offsets between adjacent continuous local 
seismic traces. So as to determine the relationships between 
the two adjacent local continuous seismic traces Sij, k (t) and 
S pq ,k(t) centered on the same sampling vertical coordinate t=t k =k 
and corresponding to different spatial coordinates i, j and p,q, 
30 there is calculated the correlation function Rij /Pq , k (h) of the two 
continuous local seismic traces. 
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The correlation function of two adjacent continuous local 
seismic traces is obtained by the following formula 



in which the numerator corresponds to the covariance 
function of Sij, k (t) and S pq , k (t) obtained by the following 



expression Cij /Pq/k (h) 

tk-A 



In this integral defining C ij/Pqfk (h), the parameter A defines an 
"investigation vertical window" about t=t k =k. 

For example, if the continuous local seismic traces Sij, k (t) 
and S P q /k (t) are trigonometric polynomials of the following form 
15 interpolating the seismic data 



3*k{t) - £a?* .cos(s<» t) + b^.sin(sai t) 
Sgjc(t) - £a^.cos(s£i?t) + bf ,sin(sfi>t) 



then it can . be shown that the covariance functions Cij, pq , k (h) is 
itself a trigonometric polynomial of the following form in which 
the. coefficients A s ij '? q ' k and - B s ij ' pq ' k depend on the coefficients 
a s ij ' k , a s pq ' k , b ij ' k and b s pq ' k : 

A conventional, and known mathematical result is that the 
correlation function 

Rit«,UW C &M*fo) translates a similarity . of 
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correlated functions when this correlation function approaches 
1. 

The study of the correlation functions of the continuous 
local seismic traces S ij/k (t) and S pq , k (t) permits defining as 
optimum offset the value hij, pq , k of h corresponding to the maximum 
of the correlation function Rij, pq ,k(h), which is to say to the 
maximum of the correlation function nearest 1. Contrary to the 
existing art based on a discrete formulation of the correlation 
function Rij, pq ,k(h), the use of a continuous formulation of the 
seismic traces and hence of the correlation function, permits 
obtaining an optimum offset which is not required to be a whole 
number multiple of the sampling interval along the vertical axis 
corresponding to the variable t. 

The definition of the optimum offset hij /Pq , k of two adjacent 
continuous local seismic traces permits obtaining a first 
approximation of the horizons passing through this reference 
continuous local seismic trace Sj.j, k (t) as shown in Figures 3 to 
5. Contrary to the prior art, the fact of using optimum offsets 
which are not required to be a whole number multiple of the 
sampling interval along the vertical axis, permits avoiding - 
errors known under the term "aliasing". 

In Figures 3 to 5, there is defined a conditional 
neighborhood Nij, k (r) of the reference local continuous seismic 
trace S±j, k (t) which is a sub-neighborhood of the initial 
neighborhood of the continuous local seismic traces taken for 
the computation of correlation and definition of the optimum 
offsets. The conditional neighborhood Nij, k (r) is selected such 
that, for any local continuous seismic trace S p£J/k (t) belonging to 
Nij, k (r), the corresponding optimum offset hij, pq/]c is such that the 
correlation R±j, pq , k (h ij/Pq/ k ) between Sij, k (t) and S pq , k (t) is greater 
than a predetermined threshold r comprised between 0 and 1. 
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As can be seen in Figure 5, the fixing of the threshold r 
also permits including in a conditional neighborhood a 
geological discontinuity or fault, which constitutes an 
important advantage relative to the prior art. 
5 The mentioned arrangements thus provide a continuous 

modeling permitting the practice of the invention. 

With reference to Figure 6, a process for smoothing or 
cleaning a subsurface in a geological structure represented by 
■seismic measurements,, comprises a first step 100 of digital 
10 modeling, as described with reference to Figures 1 to 5, to 
define the continuous local seismic traces of a multi- 
dimensional seismic matrix, calculating the optimum offsets of 
adjacent continuous local seismic traces and defining the 
conditional neighborhoods of the reference central continuous 
15 local seismic traces. 

After this first step 100 of digital modeling, a step 101 
permits selecting the property to be smoothed on a conditional 
neighborhood. This property to be smoothed on a conditional 
neighborhood can consist of any multi-dimensional function of 
20 the surface to be analyzed. 

In particular, the property to be smoothed can simply 
consist in the amplitude of the seismic signal reflected by the 
horizons of the subsurface in the direction of the measuring 
geophones . 

25 In step 102, the properties of the conditional neighborhood 

are replaced by the properties offset by the value of the 
optimum offset hij, pq , k , so as to constitute a set of properties 
relative to the same horizon. 

Thus, in the case of a reference central continuous seismic 

30 local trace Sij, k (t), all the values of S pq , k (t) of the conditional 
neighborhood are replaced by the values S pq , k ( t+hij, pq , k ) . 
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There is then carried out in step 103 an averaging of the 
offset values S pq , k (t+hij, pq , k ) established in step 102 and there is 
selected as the value of the property at the point of spatial 
coordinates i,j, and of the temporal depth coordinate t=t k =k, the 
5 average thus calculated. 

The invention is applicable to several modifications for 
establishing averages: one could for example use a weighted or 
non-weighted average; or, alternatively, one could also use an 
arithmetic average, a geometric average, a harmonic average, 

10 weighted or un-weighted . 

In the case of a weighted average, there could desirably be 
used as the weighting factor the maximum value of the 
correlation function Rij, pq ,k (h ij/Pq , k ) by thus reflecting the 
degree of confidence that can be given to each optimum offset 

15 corresponding at a maximum to a correlation function. 

The smoothing or filtering process of a property according 
to the present invention thus permits improving the contrast and 
the precision of the subsurface properties in a multi- 
dimensional seismic matrix. 

20 If necessary, so as to improve smoothing, the set of the 

processes described by the present invention can be repeated on 
the smoothing obtained during a preceding iteration. Step 104 
constitutes an iterative step in which the operator can define 
the number of iterations suitable for the resolution of his 

25 problem. 

The invention is preferably practiced on a device 
comprising means to accomplish successively the steps of the 
process described with reference to Figure 6. 

In particular, this device comprises means for visualizing 
30 seismic parameters or properties obtained in step 105 and means 
for program memorizing and intermediate computations. 

10 



Preferably, this device is a programmable device controlled 
by a computer software package, produced to practice a process 
according to the invention. 

The invention described with reference to various 
5 particular objects is not in any way thereby limited but on the 
contrary covers any modification of form or any modification of 
embodiment within the spirit and scope of the invention, the 
essential point being to use values of properties that have 
failed relative to the predetermined optimum offsets as values 
10 corresponding to the maximization of the correlation functions 
of local seismic traces that are continuous or even discrete. 
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